Magnetic field induced localization in carbon nanotubes 
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The electronic spectra of long carbon nanotubes (CNTs) can, to a very good approximation, be 
obtained using the dispersion relation of graphene with both angular and axial periodic boundary 
conditions. In short CNTs one must account for the presence of open ends, which may give rise to 
states localized at the edges. We show that when a magnetic field is applied parallel to the tube axis, 
it modifies both momentum quantization conditions, causing hitherto extended states to localize near 
the ends. This localization is gradual and initially the involved states are still conducting. Beyond 
a threshold value of the magnetic field, which depends on the nanotube chirality and length, the 
localization is complete and the transport is suppressed. 

PACS numbers: 73.63.Fg, 75.47.-m, 73.23.Ad, 85.75.-d 



The existence of geometry-induced localized states at 
the zigzag edge of graphene nanoribbons has been pre- 
dicted some years ago [1, 2], recently seen experimentally 
and shown to influence the transport in graphene quan- 
tum dots [3]. Similar states have been observed at the 
ends of a single- wall chiral nanotube studied in [4]. In 
zigzag- armchair nanotube junctions, the interface states 
calculated to appear at the junction were identified with 
the end states of the zigzag nanotube fragment [5]. 

In this Letter we predict the occurrence of localized 
states in CNTs, which is entirely due to the presence of 
a parallel magnetic field. Above a threshold flux (j)ioc 
(see Eq. (21)) these states decay exponentially with the 
distance from the nanotube end. They appear even in 
those chiral CNTs which have no localized end states 
when the magnetic field is absent. They can be found by 
a purely analytical method based on the Dirac equation 
in graphene, and the resulting energy spectrum is the 
same as that obtained by the numerical diagonalization 
of the full nanotube Hamiltonian. 

The model. Our starting point is the tight-binding 
Hamiltonian for a honeycomb lattice with one Pz orbital 
per atom and with the interatomic potential V . If we cal- 
ibrate our energy scale so that the on-site energies vanish, 
the Hamiltonian is given by 

H = Y.U,\zj){zi\, (1) 

where i and j are the lattice site indices, \zj) is a Pz 
orbital at site j and tij = {zj\V\zi) is the hopping inte- 
gral between the sites. This Hamiltonian nicely captures 
the properties of flat graphene and CNTs. In order to 
properly describe finite size nanotubes in magnetic field 
it is necessary to include the Peierls phase and curvature 
effects in the hopping elements tij. We follow here the 
approach of Ando [6]. For the sake of clarity we shall 
initially neglect the spin-orbit coupling and the Zeeman 
effect, as they do not change our main conclusion. The 
spin-dependent effects will be addressed later. 
The graphene coordinate system and the relevant real 



space vectors are shown in Fig. 1(a), while the graphene 
Brillouin zone with K and K' points is shown in Fig. 
1(b). In order to find the appropriate boundary condi- 
tions and eigenstates of CNTs we use an approach based 
on the Dirac equation treatment [2, 7]. 




FIG. 1: (a) Fragment of a graphene lattice. When consid- 
ering a CNT with chiral angle ^, we shall use the system of 
coordinates defined by directions perpendicular {x±_) and par- 
allel (x||) to the tube axis. (6) Brillouin zone of graphene. (c) 
Real and imaginary solutions of Eq. (17), determining the 
quantization of Atji as a function of n'^. (d) Some of the real 
and imaginary solutions of Eq. (17) for the Fermi subband 
{k^ = 0) of an (18,0) CNT with 100 unit cells, with and 
At 11 as functions of the magnetic flux (j). 

Parallel magnetic field. The magnetic field modifies all 
hopping integrals by a Peierls phase factor. Its form can 
be derived using the substitution p ^ p — eA and reads 

t.jiB) = t.jiO) exp 1^ ^ ' A(r) • drj . (2) 

In the cylindrical coordinates (r, (/:?,z), with the z direc- 
tion aligned with the axis of the nanotube, a parallel 
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magnetic field has coordinates B = (0, 0, 5). In the 
tangential gauge this gives A = (0, 5r/2, 0). The Peierls 
phase then becomes 



ie n 

^ Jr. 



A(r) • dr 



(3) 



where (j) is the magnetic flux threading the nanotube, 
<po = h/e the flux quantum, and cpj — cpi is the difference 
between the angular coordinates of site j and site i. 

Curvature. In a nanotube the a bonds are not orthog- 
onal to the Pz orbitals and the hopping integral tij can 
be expressed as 

tij{0) = {zj\V\zi) = K n^n • n_^-n + • Ujt, (4) 

where V-j^ and Va are hopping parameters for the corre- 
sponding bonds [6]. In our calculations we shall use the 
parameters from [8], Va = 6.38 eV and K- = —2.66 eV. 
The vector is a unit vector normal to the nanotube 
surface at the site i. The components n^n (normal) and 
n^t (tangential) are defined with respect to a plane con- 
taining the a bond between i and j and parallel to the 
CNT axis. The hopping integral tij{0) then reads 

{zj\V\zi) = cos{(pi - (pj) 
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- {Va-V^)^[l-CO^{ipi-ifj)\ , (5) 

where R is the nanotube radius and ac = 1.42 A is the 
bond length in graphene. 

In order to find the CNT spectrum it is convenient to 
express the Hamiltonian in the Bloch wave basis. The 
Bloch waves for the CNT sublattice p are given by [9] 



l^.(k)) 
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(6) 



where N is the number of the unit cells. The Bloch 
wave on the whole lattice is a linear combination of Bloch 
waves on individual sublattices and can be written as 



|$(k))= r,,(k)|$p(k)) 



In this basis the Hamiltonian acquires the form 




(7) 



^(k) = 







Hab{\^) \ 
J 



(8) 



where i^^B(k) = Xl^^i e*^"^% are the vectors con- 
necting an A sublattice atom with its neighbours, as 
shown in Fig. 1(a), and t^'s are the hopping integrals 
between an atom on sublattice A and its neighbours. Be- 
cause both the magnetic field and the curvature are uni- 
form along the whole nanotube, the t^'s do not depend 



on the position of the initial A atom. This Hamiltonian 
can be further expanded around the K {r = 1) and K' 
(r = — 1) points (see Fig. 1(b)), yielding 
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e^^^(r<+mj|) 




(9) 

where is the CNT chiral angle, indices (±, ||) denote the 
components perpendicular and parallel to the nanotube 
axis respectively, hvp = 3 |K-| <^c/2, = k — rK and 
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R ^o' 



(10a) 
(10b) 

The last term in (10a) is the Aharonov-Bohm contribu- 
tion while Akj_ and A/cy are due to the curvature, 

AA^i = ^ (l + I^V^) (11a) 
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sin(36'). (lib) 



In this derivation we used a small angle approximation, 
sm{ipi — (fij) ipi — ipj, which is good already for CNTs 
with R > 5A. The energy eigenvalues of the Hamiltonian 
are 



E± = ±hvF + (k||)2, E± := E±/{hvF). 

(12) 

Eigenf unctions of the Hamiltonian. The energy eigen- 
states are a linear combination of Bloch waves. Since 
we have expanded the Hamiltonian around the K and 
points, the corresponding Bloch waves and the coef- 
ficients r]p acquire the index r. We shall be using 

^rp{r,^) = {r\%{TK^^)). (13) 

Angular boundary condition. The wave function in the 
angular direction must be periodic. This imposes 



R' 



(14) 



which is the standard quantization condition [9]. 

Axial boundary condition. The wave function at the 
ends of the nanotube must satisfy open boundary con- 
ditions. We shall derive them for a zigzag nanotube 
{0 = 0°), but they are valid for any other chirality except 
armchair [10], provided the nanotube edge is a so-called 
minimal boundary (there are no atoms with only one 
neighbour) . 

The Hamiltonian (9) acting on the wave functions 
\^r{f^))^ (7) with (13), gives two equations: 
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We choose then, up to a normahzation factor, rjrA{i^) = 
r/^:^ + mjl and r]rB{i^) = ^IVta] e"*'^^. We can see 
from (10b) and (12) that the energies of states with /^y 
and —{i^w + 2A/cy) are the same. The energy eigenstate 
is therefore a hnear combination of both: 



V^r(r,^±) 



ai <l>^(r, (/^x,^||)) 
a2 $r(r, (/^±, 



■2Ak?^). (15) 



From the structure of the lattice in Fig. 1(a) we see that 
when the graphene patch is rohed in order to create a 
zigzag nanotube, the lower CNT edge is formed entirely 
by B sublattice atoms while the upper edge only by A 
sublattice atoms. Therefore the wave function on this 
patch must vanish at the "missing" A atoms below the 
lower edge (xy = 0) and B atoms above the upper edge 
(x|| = L). The conditions for the sublattice components 
of '0r(r, E±) are 



V^rA((x±,0),^±) =0 



(16a) 



L-iA^L ^-iK\L-i 

These equations lead to a constraint on the values of n'u , 



ai e'^'\\^-'^^\\^ + a2 g-^^-^^^S^ = 0. (16b) 



(17) 



Thus the allowed values of n\\ depend on n'^^ and in par- 
ticular on the Aharonov-Bohm flux cj). The quantity n'^^ 
can be either real or imaginary. If it is real, the wave 
function describes an extended state. If n'^^ is imaginary, 
then n\\ must be complex, with its real part equal to the 
second term in (10b). The equation (17) has then one 
trivial {n'^^ = 0) and two non-trivial solutions. The latter 
describe evanescent waves localized near the ends of the 
nanotube, because the factor exp[i(rK + a>:) -Ypi] from (6) 
acquires a damping real part. 

The regions where n'^^ is real or imaginary are determined 
by the value of rn'^ (see Fig. 1(c)). The two localized 
state solutions exist if 



states tends to with increasing magnetic flux because 
for ^ oo, k'^ for the K point solutions (see 

Fig. 1(d)). The CNT spectrum may contain localized 
{E = 0) states even at = 0, as can be seen in Fig. 
2(a) for an (18,0) CNT. If the higher {k^ 7^ 0) subbands 
lie on the Dirac cone and the condition (18) is fulfllled, 
then the lowest n\\ states in the neighbouring subbands 
are localized, while the remaining ones have energies in 
a higher range, appropriate for their subband. Whether 
the other subbands lie on the Dirac cone depends on the 
chirality and diameter of the CNT. 
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FIG. 2: Magnetic field induced localization in a (18,0) zigzag 
CNT with 100 unit cells (L = 42.6 nm, = 0°). (a) The 
spectra close to the Fermi level obtained by a numerical di- 
agonalization of the real space Hamiltonian (1) and analyti- 
cally from the Dirac- like dispersion (12) with = (^ 7^ 
states) and — [E = states), where is defined 

by (17), neglecting the spin. The black dashed line marks 
the onset of the localization, {b) The setup used for the con- 
ductance calculation, (c) Spectra obtained analytically with 
the electron spin included through (19). {d) Greyscale plot of 
conductance, in units of conductance quantum Gq — 2e^//i, 
as a function of (j) and the chemical potential E^ including 
spin effects. 



for K : hi'^> 1/L, for K' : k.'^ <-l/L. (18) 

The spectrum of the CNT is then determined by the value 
of the magnetic fleld, which enters into n'^ via (10a). In 
order to calculate the energy levels, the allowed values of 
must be found from (17) for each value of n'^ sepa- 
rately, as shown in Fig 1(d) for an (18,0) zigzag CNT. 
The analytical method described above gives a remark- 
able agreement with the spectra obtained by the numer- 
ical diagonalization of the nanotube Hamiltonian (1), 
see Figs. 2(a) and 3(a). The energy of the decaying 



Spin effects. With spin, the Bloch waves (7) become 
4-component spinors and both the spin-orbit coupling 
(SOC) and the Zeeman effect must be considered. They 
will be treated in detail elsewhere. Here we just note that 
SOC can be taken into account by yet another shift of 
while the Zeeman effect splits the energy: 

^± ^ /^l+aAA:5o, E± ^ E^ ^ aj^B^. (19) 
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FIG. 3: Localization in a (12,9) chiral CNT with 16 unit cells 
(L = 41.5 nm, = 25°). (a) Comparison of numerical and 
analytical spectra, neglecting the spin. In this case there are 
no localized states at = 0. The dashed line marks the onset 
of the localization, (b) Analytical spectrum including spin 
effects, (c) The amplitude of the highest valence eigenstate 
(obtained numerically) at each atom, projected onto x\\, for 
different values of and neglecting the spin, (d) Conductance 
as a function of magnetic flux and chemical potential 
including spin effects. 

where a = +1/ — 1 for spin parallel/antiparallel to the 
CNT axis, /i^ is the Bohr magneton and (5 is a parameter 
defining the SOC strength. In our calculations we take 
5 ~ 2.8 • 10~^, as e.g. in [11]. The resulting spectra of a 
(18,0) and (12,9) CNTs are shown in Figs. 2(c) and 3(b). 

Localization. Equations (18) define the localization 
flux 0^oc, at which the extended solutions morph into 
localized states. This threshold flux depends on the spin 
via (19) and using it together with (10a) we obtain 

(t^ioc ^^^(^J^~ ^^1^ 00, ^foc = ^loc + crRAkso 00. 

(21) 

The value of (j)ioc depends on the length of the nanotube. 
For sufficiently long CNTs the spectra are very close to 
those of the infinite nanotubes. 

The localization induced by the magnetic field is grad- 
ual, in principle allowing the involved states to conduct 
as long as the two sublattice wavefunctions overlap. The 
evolution of an eigenstate in increasing magnetic field is 
shown in Fig. 3(c), through a sequence of plots of the 
wave function amplitude at each atom, projected onto 
The apparent continuity of the curves is due to the 



overlap between close plot points; near the CNT ends 
the wave function oscillates with the azimuthal angle (f 
and the individual points can be seen clearly. Initially 
(0 = lO~^0o) the state is extended; when the magnetic 
flux reaches 0^oc5 its wave function begins to be described 
by an imaginary solution of (17). With magnetic field in- 
creasing further, the wave function decays exponentially 
with the distance from the CNT ends, the localization 
becomes complete and the state ceases to conduct. 
The above analysis is confirmed by conductance calcu- 
lations, with the CNT in a setup shown in Fig. 2(b). 
We derive the elastic linear response conductance via the 
Fisher-Lee formula for the quantum mechanical trans- 
mission: G=^ TiiTLGTRg}, where Tl/r = i{^L/R- 
^l/r)^ ^l/r is the self energy of the left or right lead 
respectively, and Q is the Green function of the cen- 
tral region dressed by the electrodes. For simulating 
bulk metal electrodes we consider wide band leads, i.e. 
J2wb{E) = —ilinS,{EF). The results shown in Figs. 
2(d) and 3(d) were obtained with ^wb = — i0.22eV. 
In both we see a gradual drop of the conductance of the 
highest valence and lowest conduction spin states, as they 
become localized in the increasing magnetic field. The 
"native" end states of the (18,0) CNT, localized also at 
= 0, can be seen in the spectrum in Fig. 2(c), but don't 
contribute to the conductance, as we expect. The good 
matching of analytical spectra of isolated CNTs and the 
conductance peaks (compare Figs. 2(c),(d) and 3(b),(d)) 
implies that even with rather strong coupling the CNT 
is sufficiently distinct from the leads for the transport to 
be determined by the spectrum of an isolated nanotube. 
The magnetic field Bioc corresponding to (j)ioc depends 
on the nanotube length and radius. For our choice of 
and Bloc of CNTs with = 7 A and L = 40 nm 
ranges from 50 T ((9 ^ 30°) to 85 T ((9 = 0°). How- 
ever, for the same nanotubes with L = 500 nm the value 
of Bloc drops to 4-42 T. Hence the localization induced 
by the magnetic field might be detected in currently ac- 
cessible transport experiments or by STM spectroscopy 
revealing localized states at the CNT ends. Moreover, 
the localization induced by a magnetic flux appears to 
be a chirality- independent phenomenon, to which only 
armchair CNTs are immune. 
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